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ABSTRACT 

The origin of the rapid quasi-periodic variabilities observed in a number of 
accreting black hole X-ray binaries is not understood. It has been suggested that 
these variabilities are associated with diskoseismic oscillation modes of the black 
hole accretion disk. In particular, in a disk with no magnetic field, the so-called 
g-modes (inertial oscillations) can be self-trapped at the inner region of the disk 
due to general relativistic effects. Real accretion disks, however, are expected to 
be turbulent and contain appreciable magnetic fields. We show in this paper that 
even a weak magnetic field (with the magnetic energy much less than the thermal 
energy) can modify or "destroy" the self-trapping zone of disk g-modes, rendering 
their existence questionable in realistic black hole accretion disks. The so-called 
corrugation modes (c-modes) are also strongly affected when the poloidal field 
approaches equal-partition. On the other hand, acoustic oscillations (p-modes), 
which do not have vertical structure, are not affected qualitatively by the mag- 
netic field, and therefore may survive in a turbulent, magnetic disk. 

Subject headings: accretion: accretion disks — black hole physics — magnetic 
fields — X-ray: binaries — hydrodynamics — MHD 



1. Introduction 

In recent years, quasi-periodic variability has been observed from a number of Galac- 
tic compact X-ray binary systems. Of particular interest are several accreting black hole 
(BH) binaries which show pairs of quasi-periodic oscillations (QPOs) of fixed frequencies 
having ratios close to 2 : 3 (for example, GRO J1655-40 shows / = 300, 450 Hz; see Remil- 
lard & McClintock 2006 for a review). The origin of these QPOs is not understood. The 
fixed frequencies and frequency ratio led to the suggestion that these QPOs involve certain 
nonlinear resonant phenomenon in the disk (e.g., coupling between the radial and vertical 
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epicyclic oscillations of the disk fluid element; Kluzniak & Abramowicz 2002), but so far no 
fluid dynamical model producing these resonances has been developed (see Rebusco 2008 
and references therein). Alternatively, it has been suggested that these QPOs may arise 
from acoustic oscillations in an accretion torus (Rezzolla et al. 2003; Lee, Abramowicz & 
Kluziniak 2004), with the oscillation frequencies determined by the (unknown) radial extent 
of the torus. Perhaps the theoretically most developed model for QPOs is the relativistic 
diskoseismic oscillation model (Kato & Fukue 1980; Okazaki et al. 1987; Nowak & Wagoner 
1991; see Wagoner 1999; Kato 2001 for reviews), in which general relativistic (GR) effect pro- 
duces trapped g-mode (also called inertial mode or inert ial-gravity mode) oscillations in the 
inner region of the disk. Various theoretical issues related to this model have been studied, 
such as the role of corotational wave absorption (Kato 2003; Li, Goodman & Narayan 2003; 
Silbergleit & Wagoner 2007) and resonant mode excitations due to global disk deformation 
(Kato 2008; Ferreira & Ogilvie 2008). 

The studies of the oscillation modes of disks/torii, such as those mentioned above, usu- 
ally assume that the unperturbed flow is laminar and has no magnetic field. Real accretion 
disks, on the other hand, are highly turbulent due to the nonlinear development of mag- 
netorotational instability (MRI) (see Balbus & Hawley 1998 for a review). The question 
therefore arises as to how the MRI-driven turbulence affects the oscillation modes obtained 
from hydrodynamical models and to what extent these trapped modes remain "valid" in a 
realistic situation. Arras, Blaes & Turner (2006) attempted to address this issue by carrying 
out MHD simulations in the shearing-box geometry. They showed that axisymmetric stand- 
ing sound waves give rise to distinct peaks in the temporal power spectrum, while inertial 
waves do not. The discrete frequencies obtained by them were due to the imposed periodic 
boundary conditions adopted in the simulations, and not due to any relativistic effect. Arras 
et al. suggested that their result posses a serious problem for QPO models based on g- 
modes (inertial waves). Recently, Reynolds & Miller (2008) reported on the results of global 
simulations of BH accretion disks (using Paczynski-Wiita pseudo-Newtonian potential) and 
showed that, while axisymmetric g-mode oscillations manifest in the hydrodynamic disk with 
no magnetic field, they disappear in the magnetic disk where MHD turbulence develops. 

In this paper, we study analytically the effects of magnetic fields on the relativistic 
diskoseismic modes in accretion disks around BHs. We consider both poloidal and toroidal 
fields and use local analysis of the full MHD equations to examine how the magnetic field 
changes the radial wave propagation diagrams for various modes. We show that the trapping 
region of g-modes can be easily "destroyed" even when the disk field strength is such that 
the associated Alfven speed is much smaller than the sound speed. On the other hand, 
the propagation characteristics of p-modes (acoustic oscillations) and c-modes are largely 
unchanged. Note that since we assume laminar flows for our unperturbed disks, we do not 
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directly address the effects of turbulence on disk modes. However, we believe that our work 
is relevant to this issue, since magnetic fields naturally arise in a turbulent disk. 

We summarize the basic MHD equations in §2 and review the properties of diskoseismic 
modes important for our analysis in §3. We examine in §4 and §5 the effects of poloidal field 
and toroidal field on those modes, respectively, and discuss the implications of our result in 
§6. 



2. Basic Equations 

We consider a non-self-gravitating accretion disk, satisfying the usual ideal MHD equa- 
tions: 

^ + V-(pv) = 0, (1) 
d 11 

-^- + (v-v)v = — vn-v$ + -r, (2) 

at p p 

— = Vx(vxB). (3) 
Here p, P, v are the fluid density, pressure and velocity, $ is the gravitational potential, and 

are the total pressure and the magnetic tension, respectively. The magnetic field B also 
satisfies the equation V • B — 0. We assume that the fluid obeys the barotropic equation of 
state P = P(p). 

We adopt the cylindrical coordinates (r, 0, z) which are centered on the central BH 
and have the z-axis in the direction perpendicular to the disk plane. The unperturbed 
background flow is assumed to be axisymmetric with a velocity field v = rfi(r)0, and 
magnetic field B = 5^,(r)0 + B z z, i.e., B z is constant while has a radial dependance. 
Force balance in the unperturbed flow implies 

G = -VU--T = Q 2 rf- V$. (5) 
P P 

Consider perturbations of the form e tm< i ) - wt . The linearized fluid equations are 

1 o %vnp o 

-iuj5p + - — (pr5v r ) H 8v^ + —{pSv z ) = 0, (6) 

r or r Oz 
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-iu5v r - = G r — - --^-5U + -(8T) r , (7) 

p par p 

-iuSv^ + ^-5v r = -—511 + -{5T) h (8) 
2\l pr p 

-iuj5v z = G z — - ~m + -(ST) X , (9) 
p poz p 

-iu5B r = ( + B g -P\ Sv r , (10) 



\ r dz / 

d d d dQ, 

-iuSB^ = - — (B^Svr) + Bz—dvj, - B 4> —5v z + r—5B r , (11) 

-en B x d imB z imB^ 

-iuooB z = —lrdv r ) ova, H -ov z , (12) 

r or r r 



where 



uj = uj — m£l (13) 

is the comoving wave frequency, and 

r dr 

is the radial epicyclic frequency. 5p, 5U, Sv, <5B are Eulerian perturbations, and p, B refer 
to the unperturbed flow variables. In addition, for barotropic fluid, we have 

bp = \5P = \(5U - ^-B ■ SB), (15) 

where c s is the sound speed. 

To perform local (WKB) analysis, we consider perturbations with spatial dependence 
e ik r r+ik z z _ j n i eac lj n g_ orc i er approximation, we keep only the radial gradient of Q(r) and 
B<f,{r) while assuming that the variation scales of all the other background quantities are 
much larger than the wavelength of the perturbation, i.e., k r , k z ^> 1/r. The linearized 
MHD equations then reduce to 



^511 + ik r 5v r + ik&Svs + ik z 5v z + «5B& + ~5B Z = 0, (16) 

pc l s 4i\pc l s 4ixpc l s 

-—511 + iuSvr + 2QSv, + ( + l h^±] SBr - ^-SBs = 0, (17) 
p \ 4np 4ixp J 2npr 

-**m - g SVr + iaSv , + (lp*t SB , + {'Mi + ■>M±) SBt = o, (18) 

p 2\l An pr \ An p 4np J 

-^5U + iu,6v, + (*M± + l M±] SBz = 0, (19) 
p \ 4np 4np J v ' 
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[ik^B^ + ik z B z )5v r + iu5B r = 0, (20) 
ikrB^Svr — ikzBzdv^ + ikzB^Svz — pQ5B r — iuSB^, = 0, (21) 
ik r B z Sv r + ik^BzSv^ — ik^B^bv z — iu5B z = 0, (22) 

where k$ = m/r. We have assumed B^ ~ r q and p = <ilnf2/(imr. Note that in deriving 
eqs. (16)-(22), we have dropped the terms proportional to G r and G z in eqs. (7) and (9): 
since G r = (f2 2 — Q 2 K )r ~ il 2 K r(H/r) 2 and G z ~ Q? K z (where Vl k is the Keplerian frequency, 
i.e., the angular frequency in the absence of pressure force), G r is much smaller than the 
other terms in eq. (7) provided that k r r 3> 1 + v^/c 2 , and G z is also negligible if we focus 
on the mid-plane of the disk. 



3. Hydrodynamic Limit: Diskoseismic Modes 

In the absence of magnetic fields, for k z , k r ^> k^, the perturbed MHD equations (16)- 
(22) lead to the dispersion relation: 

(u 2 -K 2 )(u 2 -k 2 z cl) = k 2 c^ 2 . (23) 

For k z = (or Co 2 3> k 2 c 2 ), this becomes Co 2 = k 2 c 2 s + k 2 , the usual dispersion relation for 
spiral density wave; for uj 2 k 2 c 2 , this becomes uj = ±nk z /(k 2 + k 2 ) 1 ^ 2 , describing inertial 
oscillations (e.g., Goodman 1993). 

For an accretion disk, with scale height H r, the vertical dependence of the pertur- 
bation is not well described by the plane wave e lkzZ unless k z H 3> 1. Okazaki et al. (1987) 
showed that for a thin disk, the perturbation equations can be approximately separated in r 
and z (see also Nowak & Wagoner 1991, 1992; Ipser 1994). For example, for vertically isother- 
mal disks with constant scale height H, one finds 5P(r, z), 5v r (r, z), 5v<f,(r, z) oc H n (z/H), 
while 5v z {r, z) oc H' n (z/ H), where H n (Z) (with n — 0, 1, • • •) is the Hermite polynomials and 
H' n (Z) = dH n (Z)/dZ . With this separation of variables, Okazaki et al. (1987) obtained the 
dispersion relation for a given n: 

(u 2 -K 2 )(u 2 -nnl) = k 2 c 2 u 2 , (24) 

where Q± is the vertical epicyclic frequency and is related to H by H = c s /Qj_. An important 
property of relativistic disks around BHs is that k is non-monotonic. Three types of trapped 
modes can be identified (see Fig. 1; see also Wagoner 1999, Kato 2001 and Ortega-Rodriguez 
et al. 2006 for reviews): 

(i) P-modes. For n = 0, waves can propagate in the region where u 2 > k 2 . These 
are acoustic waves (modified by disk rotation), and have also been termed inertial- acoustic 
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modes. If waves can be reflected at the disk inner radius (risccb the inner-most stable 
circular orbit), discrete p- modes can be self-trapped at the inner-most region of the disk (see 
Fig. la, lc). 

(ii) G-modes. For n > 1, waves can propagate in the region where u 2 < k 2 < nVt\ or 
Co 2 > Q 2 j_ > k 2 (note that k < Q± in GR). The former specifies the g-mode propagation zone: 
self-trapped g-modes can be maintained in the region where k peaks (for m = 0: see Fig. lb) 
or in the region where Q — n/m < uo/m < Q + n/m (for m^O: see Fig. Id) fl Because these 
discrete, self-trapped modes do not require special boundary conditions (e.g., wave reflection 
at r = risco), they have been the focus of most studies of relativistic diskoseismology. Note 
that although we call these g-modes (following the terminology of Kato 2001 and Wagoner 
1999), they have no relation to gravity waves, which are driven by buoyancy. Instead, these 
modes describe inertial oscillations, and have also been termed inertial modes (or inertial- 
gravity modes). 



1 Note that non-axisymmetric g-modes with uj/m < ri(nsco) contain corotation resonance in the wave 
zone, leading to strong damping of the mode (Kato 2003; Li et al. 2003; Zhang & Lai 2006). On the other 
hand, modes with fi(nsco) < w/m < max(f2 + n/m) do not suffer corotational damping, and are therefore 
of great interest. 
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Fig. 1. — Wave propagation diagram showing various trapped modes in BH accretion disks: 
(a) axisymmetric p-mode; (b) axisymmetric g-mode; (c)non-axisymmetric p-mode; (d) non- 
axisymmetric g-mode and c-mode. The curves depict various critical frequencies (k, y/nQ±, 
Q, Q ± k/tu, Q ± ^/nVt^/m), the vertical dotted lines denote the inner-most stable circular 
orbit (ISCO). The curvy horizontal lines specify wave propagation zones and the height of 
the line is to (for the top panels) or uj/m (for the bottom panels) of the mode. 

(iii) C-modes. For n > 1 and m > 1, the wave propagation condition uj 2 > nQ 2 ^ > k 2 
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leads to an additional wave trapping region, where u/m < Q — y/nQ±/m (see Fig. Id). Note 
that for spinning BHs, Q± < Q. Clearly, these modes exist only when Q — y^nVt±/m > 
and wave reflection occurs at r = risco Following the previous works (e.g., Kato 1990 and 
Silbergleit et al. 2001, who focused on the "fundamental" n — m — 1 mode, corresponding 
to the Lense-Thirring precession of the inner disk), we call these ("corrugation") c-modes. 

Comparing eqs. (23) and (24), we see that we can obtain the radial dispersion relation 
of different modes by adopting the vertical "quantization" condition k z = \fnjH in eq. (23), 
with k z = specifying p-modes. In a generic disk (e.g., when the disk is not isothermal 
vertically), the same "quantization" condition would not hold, but we still expect k z ~ 
1/H ~ fij_/c s for the (vertically) lowest-order g-mode or c-mode. In the next sections, we 
will adopt k z = ,/rjQ±/c s , with rj of order unity, when we study how magnetic fields modify 
low-order g-modes and c-modes. 

Our approach in this paper is based on Newtonian theory. GR effect can be incorporated 
into our analysis by using the Paczynski-Witta pseudo-Newtonian potential, $ = — Mj (r — 
2M). Alternatively, we could simply replace the Newtonian fl, Q±, k by their exact general 
relativistic counterparts (e.g., Okazaki et al. 1987): 



n = 



1 + a^/M/r 3 



1/2 



K = 



M{r 2 



AaM 1 ' 2 3a 2 
- 6Mr + SaM^r 1 ' 2 



3a 2 ) 



r 2( r 3/2 + a MV2)2 



1/2 



(25) 
(26) 
(27) 



(in geometric units such that G — c — 1), where a = J s /M is the spin parameter of the black 
hole. In general, > > k. In the case of a Schwarzschild BH, Q, = > k, with k peaks 
at r = 8M and becomes zero at risco — 6M. This non-monotonic behavior of the radial 
epicyclic frequency is preserved for Kerr BHs, and, as discussed above, is the key ingredient 
for the existence of trapped diskoseismic modes. 



4. Effect of Poloidal Fields 

We first consider the case of a pure poloidal field, with = 0. For k z ,k r 3> fc^, 
equations (16)-(22) then lead the dispersion relation: 

^ - [(kl + k 2 )(c 2 s + V\ z ) + k\v\ z + K 2 ]^ 
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h Z V A Z C : 



(e z + e r )(2c 



v 2 Az 



+ 



(K + K)v Az + 



dn 2 

dlnr 



dtt 2 
dlnr 

= 0, 



+ K 2 k 2 z c 2 s > 



(28) 



where iu 2 = B z j y/Anp. In the incompressible limit, this reduces to the dispersion relation 
found in, e.g., Balbus & Hawley (1991). For a given k = (k r ,k^,k z ), equation (28) admits 
three branches, corresponding to fast, slow magnetosonic waves and Alfven wave, all modified 
by differential rotation. For k z ^> k r , the Alfven branch can become unstable when k 2 v 2 Az < 
-dfl 2 /d\nr. This is the well-known MRI (e.g., Balbus & Hawley 1998). 



4.1. P-modes 



If k z = 0, equation (28) reduces to 



CO 



K 2 + k 2 r (ci + v 2 Az ). 



(29) 



This is almost the same expression as in pure hydrodynamic case {Cj 2 = k 2 + k 2 c 2 ), the 
only difference being that the sound speed is replaced by fast magnetosonic wave speed, 
\/ c s + v Az- Thus the basic property of p-modes is not affected by poloidal magnetic fields. 



for k 2 .: 



4.2. G-modes 

For a fixed k z = ^/H = y/rj£l±/c s (see §3), we can rewrite eq. (28) as an expression 

{u 2 - u 2 ){u 2 - u 2 ){u 2 - u 2 ) 



(c; + v 2 Az )k 



(u) 2 — LU 2 )(0J 2 — uj) 



The five critical frequencies are given by 



UJn 



OJn 



10 A 



I \k 2 + 2r](Q ± ) 2 b 2 + vV + I6ri(tt ± n) 2 b 2 



\2l2 



1 + b 2 ' 



k 2 + 2r](n ± ) 2 b 2 - Vk 4 + lQr](n ± n) 2 b 2 



where b = va z /c s . 



(30) 

(31) 
(32) 
(33) 
(34) 

(35) 
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Equation (30) allows us to identify various wave propagation regions (k 2 > 0). We first 
consider subthermal fields, with b < 1. When b < 0.4 (and with r] — 1 for the lowest order 
g- modes), the five critical frequencies satisfy uj\ > uj\ > > uj\ > > uj\ in the inner 
region of the disk. Thus there are three wave propagation regions: 



Region I : lj 2 > uj\, 
Region II : uj\ < lj 2 < uj\, 
Region III : Cj 2 < uj\. 



(36) 
(37) 
(38) 
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r/M 



Fig. 2. — The effect of poloidal magnetic field on the g-mode propagation zone for m = 
0, rj = 1. In each panel, the upper three curves are u>2 (eq. [32]) and the lower two curves 
are uj% (eq. [33]). Axisymmetric g- modes of frequency us can propagate in the region where 
co> 3 < uj < uj 2 . The solid line refers to the case of b = 0, the short-dashed line b = 0.03 and 
the long-dashed line b = 0.1, where b = va z /c s , with vaz — B z / ^^np (Alfven speed) and c s 
the sound speed. The vertical dotted lines correspond to the inner disk radius at ISCO. The 
curvy horizontal lines specify the wave propagation zones, and the height of the line is uj of 
the mode. The lower and upper panels are for the case of a Schwarzschild BH (a = 0) and a 
Kerr BH (a = 0.8M), respectively. The angular frequencies are in units of M _1 = c 3 / (GM). 
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Region II corresponds to the original g-mode cavity modified by the magnetic field; in 
the zero field limit, u 3 = 0, u 2 = k and eq. (37) reduces to u 2 < k 2 . Fig. 2 depicts the critical 
frequencies U2 and 003 for several values of b. This also serves as the propagation diagram 
for m = g- modes (wave can propagate in region where 003 < uo < uj 2 ). We see that as 
the magnetic field increases, the g-mode self-trapping zone gradually shrinks and disappears 
even when the magnetic field is still very subthermal (for a Schwarzschild BH, this occurs for 
b > 0.08). More precisely, the g-mode cavity can still exist for large b, but it now requires a 
reflection boundary at r IS co- This behavior can be easily understood by inspecting eq. (32): 
While n peaks at some radius r max , Qj_ and fl both increase monotonically with decreasing 
r. Since Qj_ and Q are much larger than k in the inner region of the disk, the 2(Q_ s _b) 2 term 
or the 4Q±Qb term can dominate over k 2 even when b is still small, therefore making the 
self-trapping zone disappear. Roughly, this occurs at b > 6 crit ~ (/t 2 /2f2_|_f2) rmax . 
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Fig. 3. — The effect of poloidal magnetic field on the g-mode propagation zone for m ^ 
0, rj — 1, a — 0. The three panels are for b = va z /c s = 0, 0.005 and 0.015. The solid line, 
short-dashed lines and the long-dashed lines show Q, Q ± uj 2 and Q ± u> 2 /2, respectively. In 
the bottom panel, the dotted lines show Q ± uj% (In the upper and middle panels, Q ± u 3 
almost concide with Q, since 0^3 = for 6 = and U3 <C Q for b <C 1). Non-axisymmetric 
g-modes can propagate in the region where Q — uj 2 /m < u/m < Q — 003/m or Q + u> 3 /m < 
u/m < Q + u 2 /m. The vertical dotted lines correspond to the inner disk radius at ISCO. 
The curvy horizontal lines in top and middle panels specify wave propagation zones and the 
height of the line is u/m of the mode. Note that the self-trapping zone (depicted in the 
upper and middle panels) disappears as b increases. The angular frequencies are in units of 
M~ l = c 3 /(GM). 
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For non-axisymmetric perturbations (m / 0), the wave propagation region II is deter- 
mined by (i) Q — iv 2 /m < iv/m < Q + iv 2 /m, and (ii) iv/m > Q + Lv 3 /m or uj /m < Q — iv 3 /m. 
Fig. 3 shows the propagation diagram. As mentioned before (see Footnote 1), for 6 = 0, 
only the modes with iv > mf2(risco) are of interest, since otherwise there is a corotation 
resonance in the wave zone, leading to strong mode damping (Kato 2003; Li, Narayan & 
Goodman 2003; Zhang & Lai 2006). Thus, self-trapped g- modes reside around the radius 
where Q + iv 2 /m is the maximum (and this maximum arises because k depends nonmono- 
tonically on r). We see from Fig. 3 that this g-mode self-trapping region disappears as b 
increases. The larger m is, the more fragile is the cavity. For example, the m = 1 cavity 
disappears for b > 0.015, while for m = 2, this occurs for b > 0.005. 



4.3. C-modes 



For m/0 and rj ~ 1, equation (30) also describes trapped c-modes. When b < 0.4, io\ is 
the largest among all the critical frequencies and the c-mode propagation zone corresponds 
to Region I (see eq. [36]). Note that since ui is not affected by the magnetic field, the 
trapping region is determined by u/m < Q — ui/m = Q — Q±/m (for rj — 1, see the upper 
panel of Fig. 4; cf. Fig. Id). When b > 0.4, the ordering between ui and uo 2 switches and 
c-modes propagate in the region where uj 2 > u 2 , with the trapping zone determined by 
iv /m < Q — iv 2 /m (see the bottom panel of Fig. 4). Thus, in the presence of a reflection 
boundary at risco; c-modes are not affected by the poloidal magnetic field when b < 0.4, 
but can be appreciably modified when b > 0.4. 



- 15 - 




r/M 



Fig. 4. — The effect of poloidal magnetic field on c-modes with m — 1 and a = 0.2M. The 
upper panel shows the original c-modes since b is small and uj\ = Q±; in the bottom panel, 
with a large b, the c-mode trapping zone is instead bounded by the inner reflection boudary 
and Q — The vertical dotted line refers to the inner disk radius at ISCO. The curvy 
horizontal lines specify wave propagation zones and the height of the line is u of the mode. 
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Fig. 5. — Wave propagation diagram for non-axisymmetric g- modes and c- modes, for b = 
VAz/c s = 0.7, m — 2, r) — 1, and a = 0. The solid line, dot-short dashed lines, long-dashed 
lines, dotted lines and short-dashed lines show Q, Q±uj^/2, Q±uj 3 /2, fl±a>i/2 and Q±uj 2 /2 
(eqs. [31]- [34]), respectively. The vertical dotted line shows the inner disc radius. The curvy 
horizontal lines specify wave propagation zones and the height of the line is u/2 of the mode. 
The singular points (where k r — > oo) are indicated by filled circles. The angular frequencies 
are in units of M -1 = c 3 /(GM). 



-17- 



From eq. (30) we can identify other wave propagation zones (see eq. [38]). Fig. 5 
gives an example, for m = 2, b = 0.7. Note that, except for the c-mode trapping zone 
discussed above, all the propagation zones are bounded by at least one "singular point" 
(where k r — > oo). Unlike the turning point (k r — ► 0) associated with wave reflection, wave 
absorption is expected to occur at these singular points (see Zhang & Lai 2006; Tsang & 
Lai 2008a and references therein). Thus, the new wave trapping regions given by eq. (38) 
will not lead to interesting global oscillation modes (Note that in the case of b = 0.7, the 
ordering of five critical frequencies is different from the one described in §4.2. However, our 
conclusion still holds true, i.e., there is no chance to form a wave trapping zone bounded 
by two reflection points other than the c-mode oscillation region, which is bounded by a 
reflection point and the ISCO). 



5. Effect of Toroidal Fields 

In this section, we consider the effect of a pure toroidal field, with B z = 0. Various 
instabilities may exist for such field geometry, depending upon the rotation profile Q(r) and 
the magnetic field profile B^(r) (e.g., Acheson & Gibbons 1978; Terquem & Papaloizou 
1996). Here we focus on how affects the diskoseismic modes. 



5.1. P-modes 

With k z = 0, equations (16)- (22) reduces to 

& = * 2 + + v%), (39) 

where VA<j> = B^/^Airp. Thus, the toroidal field affects p-modes in the same way as the 
poloidal field does (see §4.1). 



5.2. G-modes 



Since the general dispersion relation for m ^ is quite complicated, here we focus on 
axisymmetric perturbations^] With m — 0, Equations (16)-(22) lead to 



)2 Acj> c 



+2(1 - q)v 2 M c 2 s k 2 Jr 2 = 0, (40) 



2 Since c- modes necessarily require m > 0, our analysis here cannot be applied to c- modes. 
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where q = dlnB^/dlnr. Solving for k 2 , 



we have 




(uj 2 - u 2 + )(u 2 - cjj 

(C 2 s + V 2 A ^ 2 



(41) 




l0«2_^2 (1 + ^ )]2 _ 



8(1 - gVW, 



(42) 



where b<p = VA<f,/c s and we have used k z = y/rj/H = y/rjQ±/c s as in §4. Clearly, for = 0, 
eq. (42) reduces to eq. (24). 

When q — 1 (i.e., oc r), eq. (42) gives = r](fi^) 2 (l + ti^), and cj 2 = k 2 . Since a; 2 , 
is independent of B^, the g-mode propagation zone is unaffected no matter how strong the 
field is. When q ^ 1, as long as VA<f, <C (which is valid in most disk situations), the 
8(1 — q)r}v A(j) Q, 2 ± / r 2 term in eq. (42) represents only a small correction, i.e., u; 2 is still very 
close to k 2 . Thus for general toroidal field satisfying va$ ^±r, the axisymmetric g-mode 
propagation zone is not affected by the magnetic field. 



In this paper we have studied the effects of both poloidal and toroidal magnetic fields on 
the diskoseismic modes in BH accretion disks. Previous works by Kato, Wagoner and others 
have been based on hydrodynamic disks with no magnetic field. The key finding of our 
paper is that the g-mode self-trapping zone (which arises from GR effect) disappears when 
the disk contains even a small poloidal magnetic field, corresponding to va z /c s = 0.01 — 0.1 
(see Fig. 2-3; va z is the Alfven speed and c s is the sound speed). It is well-known that 
the combination of a weak poloridal field and differential rotation gives rise to MRI, making 
real astrophysical disks turbulent. Earlier numerical simulations indicated that the magnetic 
field grows as MRI develops, until it saturates at va z /c s ~ 0.1 — 1, with the toroidal field 
stronger than the poloidal field by a factor of a few (see, e.g., Hawley et al. 1996; Balbus 
& Hawley 1998). Recent simulations showed that the turbulent state depends strongly on 
the net magnetic flux through the disk (e.g., Fromang & Papaloizou 2007; Simon, Hawley & 
Beckwith 2008). In any case, it is likely that the magnetic field in a turbulent disk is large 
enough to "destroy" the g-mode self-trapping zone. 

Thus, the g-mode properties (including the frequencies and excitations) derived from 
hydrodynamical models are unlikely to be applicable to real BH accretion disks. The dis- 
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appearance of the g-mode trapping zone might also explain why Arras et al. (2006) and 
Reynolds & Miller (2008) did not see any global g-modes in their MHD simulations. 

As mentioned in §1, g-mode oscillations have been considered a promising candidate to 
explain QPOs in BH X-ray binaries. Theoretically, these modes are appealing because in 
hydrodynamic disks their existence depends on general relativistic effect and does not require 
special disk boundary conditions. Our analytical results presented in this paper, together 
with recent numerical simulations (Arras et al. 2006; Reynolds & Miller 2008), suggest that 
magnetic fields and turbulence associated with real accretion disks can change this picture 
significantly. 

While g-modes can be easily modified or "destroyed" by magnetic fields, our analysis 
showed that p-modes are not affected qualitatively. The magnetic field simply changes the 
sound speed to the fast magnetosonic wave speed and leaves the p-mode propagation diagram 
unchanged. We also showed that a weak poloidal field {vaz/c s <C 1) does not affect the c- 
mode propagation zone, although a stronger field modifies it. Our results therefore suggests 
that global p-mode oscillation is robust and may exist in real BH accretion disks, provided 
that partial wave reflection at the disk inner edge can be achievedjfl Of particular interest 
is the non-axisymmetric p-modes, since they may be excited by instabilites associated with 
corotation resonance (Tsang & Lai 2008a, b). 

We thank David Tsang for useful discussion. This work has been supported in part by 
NASA Grant NNX07AG81G and by NSF grant AST 0707628. 
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